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1 Introduction 



We have been living with QCD spectral sum rules (QSSR) (or QCD sum rules, or ITEP 
sum rules, or hadronic sum rules...) for 15 years, within the impressive ability of the 
method for describing the complex phenomena of hadronic physics with the few universal 
"fundamental" parameters of the QCD Lagrangian (QCD coupling a^, quark masses and 
vacuum condensates built from the quarks and/or gluon fields, which parametrize the non- 
perturbative phenomena). The approach might be very close to the lattice calculations as 
it also uses the first principles of QCD, but unlike the case of the lattice, which is based on 
sophisticated numerical simulations, QSSR is quite simple as it is a semi-analytic approach 
based on a semiperturbative expansion and Feynmann graph techniques implemented 
in an Operator Product Expansion (OPE), where the condensates contribute as higher 
dimension operators in the OPE. In this approach, one can really control and in some 
sense understand the origin of the numbers obtained from the analysis. 

Since its early days and up to now, the field has, unfortunately, suffered from some 
apparent controversy and (worse) from some irrationnal emotional fights between different 
groups, which obscure its beauty. But, when you start d'y fourrer votre nez , you realize 
that the QCD spectral sum rule is one of the most interesting discoveries of the last decade 
since, with its simplicity, it can describe in an elegant way the complexity of the hadron 
phenomena, without waiting for a complete understanding of the confinement problem. 

One can fairly say that QCD spectral sum rules already started, before QCD, in the 
time of current algebra in 60, when peoples proposed different ad hoc superconvergence 
sum rules, especially the Weinberg and Das-Mathur-Okubo sum rules, which came under 
control only within the advent of QCD [|l|. However, the main flow comes from the classic 
paper of Shifman-Vainshtein-Zakharov (hereafter referred to as SVZ), which goes 
beyond the naive perturbation theory thanks to the inclusion of the vacuum condensate 
effects in the OPE. More details and more complete discussions of QSSR and its various 
applications to hadron physics can be found, for instance, in f^. 

In this talk, I shall present aspects of QSSR in the analysis of the properties of heavy 
flavours. As I am limited in space-time (an extended and updated version of this talk 
will be published in the Proceedings of QCD 94 at Montpellier), I cannot cover in detail 
all QSSR applications to the heavy-quark physics here. I will only shortly discuss the 
following topics, which I think are important in the development of the understanding of 
the heavy-quark properties. These concern the determination of the: 

- heavy-quark "perturbative" pole and running masses, 

- decay constants of the D, B and D^, Bs mesons and the iJ^-parameter, 

- semileptonic and radiative decay form factors of the 5, 

- slope of the Isgur-Wise function and the emerged value of V^t, 

- properties of hybrids and Be mesons. 
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2 QCD spectral sum rules 



In order to illustrate the QSSR method in a pedagogical way, let us consider the two-point 
correlator: 

^^(g^M,2) ^ ^|d'xe^'^^ {0\TJi:ix)iJao))^\0) (1) 



where Jf^{x) = lrf'^b{x) is the local vector current of the 6— quark. The correlator obeys 
the well-known Kallen-Lehmann dispersion relation: 

r-oo fjj^ 

nb(q'^,M?)= ImnJt) + subtractions, (2) 

J4Af2 t — q^ — ie IT 

which expresses in a clear way the duality between the spectral function Im 11;, (t) , which 
can be measured experimentally, as here it is related to the e'^e~ into T-like states total 
cross-section, while Ilb{q^, M^) can be calculated directly in QCD, even at = 0, thanks 
to the fact that - > A2. The QSSR is an improvement on the previous dispersion 
relation. 

In the QCD side, such an improvement is achieved by adding to the usual perturbative 
expression of the correlator, the non-perturbative contributions as parametrized by the 
vacuum condensates of higher and higher dimensions in the OPE |Q : 

^,(g^M,2)^ E ^ E C7(^H?^M^^)(0(M)), (3) 

D=0,2,4,... (jWfc — q j dimO=D 

where /i is an arbitrary scale that separates the long- and short- distance dynamics; C'-"'^ 
are the Wilson coefficients calculable in perturbative QCD by means of Feynman diagrams 
techniques: D = corresponds to the case of naive perturbative contribution; (O) are 
the non-perturbative condensates built from the quarks or/and gluon fields. For D = 4, 
the condensates that can be formed are the quark Mi{ijjil>) and gluon {agC^) ones; for 
D = 5, one can have the mixed quark-gluon condensate {ipa ^v^"" / '^G^^ if)) , while for D = Q 
one has, for instance, the triple gluon gfabc{G"'G''G'^) and the four-quark as{ipTiipijjT2ip) , 
where Fj are generic notations for any Dirac and colour matrices. The validity of this 
expansion has been understood formally using renormalon techniques |^ and by building 
renormalization invariant combinations of the condensates (Appendix of 0]). The SVZ- 
expansion is phenomenologically confirmed from the unexpected accurate determination 
of the QCD coupling from semi-inclusive tau decays @, . In the present case of heavy- 
heavy correlator the OPE is much simpler, as one can show 0, that the heavy-quark 
condensate effects can be included into those of the gluon condensates, such that up to 
D < 6, only the G^ and G^ condensates appear in the OPE. Indeed, SVZ has, originally, 
exploited this feature for their first estimate of the gluon condensate value. 

For the phenomenological side, the improvement comes from the uses of either finite 
number of derivatives and finite values of q^ (moment sum rules): 



r°° dt 1 ^ , , 
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or infinite number of derivatives and infinite values of g^, but keeping tlieir ratio fixed as 
T = n/q^ (Laplace or Borel or exponential sum rules): 

poo 1 

C{t,MI)= (it exp(-tr) - Imnfe(t). (5) 

There also exists non-relativistic versions of these two sum rules, which are convenient 
quantities to work with in the large- quark- mass limit. In these cases, one introduces 
non-relativistic variables E and t^'- 

t = {E + Mbf and tn = ^MbT. (6) 

In the previous sum rules, the gain comes from the weight factors, which enhance the 
contribution of the lowest ground state meson to the spectral integral. This fact makes 
the simple duality ansatz parametrization: 

"one resonance"5(t - M^) + "QCD continuum"e(t - tj, (7) 

of the spectral function to give a very good description of the spectral integral, where the 
resonance enters via its coupling to the quark current. In the case of the T, this coupling 
can be defined as: 

(0|67'^6|T) = V2^. (8) 

The previous feature has been tested in the light-quark channel from the e"'"e^ ^1=1 
hadron data and in the heavy-quark ones from the e^e" ^ T or ip data, within a good 
accuracy. To the previous sum rules, one can also add the ratios: 

7^(") = -^^ and 7^. = --log£, (9) 

and their finite energy sum rule (FESR) variants, in order to fix the squared mass of the 
ground state. In principle, the pairs {n,tc), ij^tc) are free external parameters in the 
analysis, so that the optimal result should be insensitive to their variations. Stability 
criteria, which are equivalent to the variational method, state that the optimal results 
should be obtained at the minimas or at the inflexion points in n or r, while stability in 
tc is useful to control the sensitivity of the result in the changes of tc-values. To these 
stability criteria are added constraints from local duality FESR, which correlate the de- 
value to those of the ground state mass and coupling [Q . Stability criteria have also been 
tested in models such as harmonic oscillator, where the exact and approximate solutions 



are known [|T0|. The most conservative optimization criteria which include various types 
of optimizations in the literature are the following: the optimal result is obtained in the 
region, from the beginning oi r/n stability (this corresponds in most of the cases to the 
so-called plateau discussed often in the literature, but to my opinion, the interpretation 
of this nice plateau as a good sign of a good continuum model is not sufficient, in the 
sense that the flatness of the curve extends in the uninteresting high-energy region where 
the properties of the ground state are lost), until the beginning of the tc stability, where 
the value of tc corresponds to about the one fixed by FESR duality constraints. The 
earlier sum rule window introduced by SVZ, stating that the optimal result should be 
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in the region where both the non-perturbative and continuum contributions are small is 
included in the previous region. Indeed, at the stability point, we have an equilibrium 
between the continuum and non-perturbative contributions, which are both small, while 
the OPE is still convergent at this point. 



3 The "perturbative" masses 

We shall only discuss the determinations of the quark masses from relativistic sum rules 
which use a truncated perturbative series, where the perturbative masses entering in the 
sum rule analysis are then well-defined. This is not the case of non-relatisvistic sum rules 
and heavy-quark effective theory where the quark is considered to be infinitely massive 
and the QCD series is summed up. In this case, the pole mass is ill-defined, due to the 
parasitic renormalon effects, which induce after a resummation of the QCD series at large 
order a "non-perturbative" piece of the order of A (however, this effect is regular izat ion- 
scheme dependent and is absent in the MS -scheme) [0. Also, in this infinite mass 
limit, the effects of the coulombic terms are important which render the estimate from 
non-relativistic sum rules very sensitive to the value of as used in the analysis. These 
non-perturbative effects induce a kind of dressed quark where the associated ill-defined 
pole mass has a strength similar to the one used in potential models, but higher than 
the one obtained from the standard relativistic sum rule at finite values of the quark 
mass and for a truncated QCD series. In the same way, non-perturbative effects, due to 
hadronization, shift, to higher values, the mass obtained from a fit of the inclusive decays 
and deep inelastic scattering data. Within the same reasoning, we would also expect that 
the pole mass from the lattice contains a non-perturbative piece. 

The masses obtained from the relativistic sum rules should be identified, like in the case of 
the light quarks, with the so-called current mass entering in the QCD Lagrangian, which is 
a well-defined quantity in perturbation theory within the MS'-scheme. Using relativistic 
sum rules in the ip and T channels, different groups have obtained the Euclidean mass 



M,(/ = -M^) ~ (1.26 ± 0.02) GeV M6(/ = -M^) ~ (4.18 ± 0.02) GeV, (10) 

defined in the Landau gaugeQ. However, the use of this gauge- dependent mass is not 
inconvenient as we know that the correlator entering the sum rule analysis is gauge- 
invariant. Instead, the use of the Euclidean mass minimizes the size of the radiative 
corrections at the order where the sum rule analysis is done, making the result insensitive 
to the error in the value of a^. One can translate this result into the perturbative pole 
mass through the relation: 

MQip' = Ml) ^ Mq{p' = -Ml) {l + 21og2^) , (11) 



^We have included in the average the shghtly lower and more precise value Mb(jj^ — ~M^) — (4.17 ± 
0.02) GeV [12] obtained in this channel from relativistic moment sum rules. 
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from which one can deduce the value of the perturbative pole masses which are gauge- 
invariant: 

Me(p2 = M^) ~ (1.46 ± 0.05) GeV Mh{p'^ = M^) ~ (4.64 ± 0.06) GeV, (12) 

for A5 ~ (180 ± 80) MeV. One can also estimate the perturbative 6-pole mass by using 
the ratios of the relativistic sum rules in the B and B* channels. One obtains : 

Mfe(p2 = Ml) ~ (4.56 ± 0.05) GeV. (13) 

By taking the average of their values from the previous two independent sources, one can 
deduce (see also ITSf): 

Mfe(/ = Ml) ~ (4.59 ± 0.04) GeV. (14) 

The value of the perturbative 6-pole mass is definitely lower than the ones, in the range 
4.8-5.3 GeV, from non-relativistic approaches, inclusive data and lattice calculations, due 
to the presence of the non-perturbative components induced by the summation of the 
QCD series in these latter cases. The running masses in the MS'-scheme can be deduced 
directly from the Euclidian masses in (10) or equivalently from the perturbative pole 
masses in (12) and (14) through (derived for the first time in [0): 

Mq{p^ = Ml) ^ Mq{MI) (1 + . (15) 

One obtains: 

Me(l GeV) ~ (1.40 ±0.06) GeV Mb(l GeV) ~ (5.87 ± 0.06) GeV, (16) 
which, combined with the s-quark mass values in [|l^, gives: 

Mfe/m, ~ 36.7 ±2.2, (17) 
a result of a great interest for model-building and GUT-phenomenology. 



4 The pseudoscalar decay constants and Bb 

The decay constants fp of a pseudoscalar meson P are defined as: 

(m, ± MQ)(0|g(^75)Q|P) = v^M|,/p, (18) 

where in this normalisation = 93.3 MeV. A rigorous upper bound on these couplings 
can be derived from the second-lowest super convergent moment: 



2! [dq'f 



(19) 

g2=0 



where \E'5 is the two-point correlator associated to the pseudoscalar current. Using the 
positivity of the higher-state contributions to the spectral function, on can deduce [^: 



where one should not misinterpret the mass- dependence in this expression compared to 
the one expected from heavy quark symmetry. Applying this result to the D-meson, one 
obtains: 

fD < 2.UU. (21) 

Although presumably quite weak, this bound, when combined with the recent determi- 
nation to two loops [ll6| : 

^^(1.15 ±0.04)/^, (22) 

JD 



implies 



fn, < (2.46 ± 0.09)/,, 



(23) 



which is useful for a comparison with the recent measurement of fr,^ from WA75: fo^ — 
(1.76 ± 0.52)/, and from CLEO: /d, ^ (2.61 ± 0.49)/,. One cannot push, however, 
the uses of the moments to higher n-values in this D-channel, in order to minimize the 
continuum contribution to the sum rule with the aim to derive an estimate of the decay 
constant because the QCD-series will not converge at higher n-values. The estimate of 



17|, based on the lowest moment, indeed suffers from the continuum sensitivity as a little 



change in the continuum threshold makes a big change of the estimate, in such a way 
that the result becomes unreliable. In the D-channel, the most appropriate sum rule is 
the Laplace sum rule. The results from different groups are consistent for a given value 
of the c-quark mass and lead to the average ||T8|, |14| : 



/d~ (1.31 ±0.12)/, 



/z5.^ (1.51 ±0.15)/,. 



(24) 



For the 5-meson, one can either work with the Laplace, moments or their non-relativistic 
variants. Given the previous value of Mj,, these different methods give consistent values of 
fs, though the one from the non-relativistic sum rule is very inaccurate due to the huge 



effect of the radiative corrections in this method. The average value of fs is ||19|| : 



fB ^ (1.60 ± 0.26)/, 



(25) 



while 



L 

fE 



1.16 ±0.04, 



(26) 



where the most accurate estimate comes from the "relativistic" Laplace sum rule. One 
could notice, since the first result fs — fo of |TB|, a large violation of the scaling law 
expected from heavy-quark symmetry. Indeed, this is due to the large 1 /M;,-correction 
found from the HQET sum rule pp] , pT| and from the one in full QCD |T^, 



(0.42 ±0.07) GeV3/^{i_(^^^±^p^j, 



(27) 



which is due to the meson-mass gap 5M = Mb — and to the continuum energy 

[0, ^ (E, ^ f5M B): 



/ 



51 



TT 



6M 3 Er 



2M, 



+ 



{2i 
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One can notice that the apparent disagreement among different existing QSSR numerical 
results in the literature is mainly due to the different values of the quark masses used 
because the decay constants are very sensitive to that quantity as shown explicitly in ||16 



Indeed, one can also exploit this sensitivity in order to deduce the value of the quark mass 
when the experimental measurement of these decay constants will become available. 

Finally, let me also mention that we have also tested the validity of the vacuum saturation 
= 1 of the bag parameter, using a sum rule analysis of the four-quark two-point 
correlator to two loops [^. We found that the radiative corrections are quite small. 
Under some physically reasonnable assumptions for the spectral function, we found that 
the vacuum saturation estimate is only violated by about 15%, giving: 

5b ~1 ±0.15. (29) 

By combining this result with the one for fs in (25), we deduce: 

/b /b^ ~211±36Mel^, (30) 

where we have used the normalization /^r = 132 MeV, which is a/2 the one defined in (18). 
These previous results are in excellent agreement with the present lattice calculations ^4 \ . 
In particular: 

/b/b^~ 220±40 Mel/, (31) 
A tentative average of these two numbers leads to the "world average" : 

/b/b^~ 215±27 Mel/, (32) 
which is useful in the phenomenological analysis of the B^-B^-mixing. 



5 Rare and semileptonic decays 

One can extend the analysis done for the two-point correlator to the more complicated 
case of three-point function in order to study the form factors related to the B K*'j 
and — *■ p/tt semileptonic decays. In so doing, one can consider the generic three-point 
function: 

Vip,p',q') ^-Jd'x d'y e'^^'^-^y^ (0|TJi(x)O(0)4(y)|0) (33) 

where J^, Jb are the currents of the light and i?-mesons; O is the weak operator specific 
for each process (penguin for the -ft'*7, weak current for the semileptonic); q = p — p'. 
The vertex obeys the double dispersion relation : 

V{p\p'\q')c. r t ■ r ' • -,^rnV{s^s',q') + ... (34) 

J Ml s — p'^ — le Jm\ s — — le 7r^ 

As usual, the QCD part enters in the LHS of the sum rule, while the experimental 
observables can be introduced through the spectral function after the introduction of 
the intermediate states. The improvement of the dispersion relation can be done in the 
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way discussed previously for the two-point function. In the case of the heavy-to-hght 
transition, the only possible improvement whith a good Mb-behaviour at large Mf, is 
the so-called hybrid sum rule (HSR) corresponding to the uses of the moments for the 
heavy-quark channel and to the Laplace for the light one [|14|, p5| : 

1 roo POO 

n{n, r') = - 1^ ds' e-^ ^ lmV{s, s\ q'). (35) 

h 

We have studied analytically the different form factors entering in the previous processes 
p6|. They are defined as: 



2V 



Ms + Mp 

{7T{p)\u-f^b\B{p)) = f+{p + p')^^ + f-{p-p')f,, 



+ {s(M|-M^^)-e*g(p + p')4^- (36) 

We found that they are dominated, for Mb oo, by the effect of the light-quark conden- 
sate, which dictates the Mfe-behaviour of the form factors to be typically of the form: 

(3r, 

where Ip is the integral from the perturbative triangle graph, which is constant as 
t'^Ec/{dd) {t'^ and Ec are the continuum thresholds of the light and h quarks) for large 
value of Mfe. It indicates that at = 0, all form factors behave like a/M^, although, in 
most cases, the coefficient of the term is large. The study of the g^-behaviours of the 

form factors indicates that with the exception of the Ai-form factor, the g^-dependence 
of the others is only due to the non- leading (in 1/Mb) perturbative graph such that for 
Mb oo, these form factors remain constant from = to q'^ax- The resulting Mb- 
behaviour at q^ax the one expected from the heavy quark symmetry. The numerical 
effect of this g^-dependence at finite values of Mb is a polynomial in g^, which mimics quite 
well the usual pole parametrization for a pole mass of about 6-7 GeV. The situation for 
the Ai is drastically different from the other ones as here the Wilson coefficient of the 
(dd) condensate contains a g^-dependence and reads: 

which, for g^^^, = {Mb — MpY, gives the expected behaviour: 

MqLx) ~ (39) 



/Mb 
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One can notice that the g^-dependence of Ai is in complete contradiction with the pole 
behaviour, as has been also noticed in the numerical analysis of |^. It is urgent and 
important to test this feature experimentally. One can finally notice that due to the 
overall l/fs factor, all form factors have a large l/M;,-correction. 

In the numerical analysis, we obtain at = 0, the value of the B K*'-/ form factor: 

Ff"''' ~ 0.27 ± 0.03, \ - 1.14 ±0.02, (40) 

which leads to the branching ratio (4.5 ± 1.1) x 10~^, in perfect agreement with the CLEO 



data and with the estimate in One should also notice that, in this case, the coefficient 
of the correction is very large, which makes the extrapolation of the the c-quark 

results to higher values of the quark mass dangerous. This extrapolation is often done in 
some lattice calculations. 

For the semileptonic decays, QSSR give a good determination of the ratios of the form 
factors with the values p5| : 

4aM.^. 1.11 ±0.01 
Ai(0) Ai(0) 

■ ^ ' ~ 1.18 ±0.06 

(41) 

Combining these results with the "world average" value of /+(0) = 0.25 ± 0.02 and the 
one of Ff~*''{0), one can deduce the rate and polarization: 

(4.3±0.7)|K6|' X 10^2 g-i 
Tp/T^ ^ 0.9 ±0.2 

r+/r_ ~ 0.20 ±0.01 

a = 2-^ - 1 ~ -(0.60±0.01). (42) 

These results are much more precise than the ones from a direct estimate of the absolute 
values of the form factors due to the cancellation of systematic errors in the ratios. They 
indicate that, we are on the way to reach Vub with a good accuracy from the exclusive 
modes. Also here, due mainly to the non-pole behaviour of Af, the ratio between the 
widths into p and into vr is about 1, while in different pole models, it ranges from 3 to 10. 
For the asymmetry, one obtains a large negative value of a, contrary to the case of the 
pole models. 



6 Slope of the Isgur— Wise function and Vcb 



Let me now discuss the slope of the Isgur- Wise function. Taron-de Rafael |^has ex- 
ploited the analyticity of the elastic 6-number form factor F defined as: 

{Bip') |67'^6|S(6)) = {p + p'YF{q^) (43) 
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which is normahzed as F{0) = 1 in the large mass hmit Mb — M^,. Using the positivity 
of the vector spectral function and a mapping in order to get a bound on the slope of F 
outside the physical cut, they obtained a rigorous but weak bound: 

F'{vv' = 1) > -6. (44) 

Including the effects of the T states below BB thresholds by assuming that the TBB 
couplings are of the order of 1, the bound becomes stronger: 

F'{vv' = 1) > -1.5. (45) 

Using QSSR, we can estimate the part of these couplings entering in the elastic form 



factor. We obtain the value of their sum 30 



5^^^55-0.34 ±0.02. (46) 

In order to be conservative, we have considered the previous estimate within a factor 3 
larger. We thus obtained the improved bound 

F'{vv' = 1) > -1.34, (47) 

but the gain is not much, compared with the previous one. Using the relation of the form 
factor with the slope of the Isgur-Wise function, which differs by -16/75 log as(Mb) 
one can deduce the final bound: 

C'(l) > -1.04. (48) 

However, one can also use the QSSR expression of the Isgur-Wise function from vertex 
sum rules in order to extract the slope analytically. The physical IW function reads: 

Cpkysiy ^ W) = (^^^ 1 1 + ^/(^) _ {ddyg{y) + (49) 

{asG'yh{y)+g{dGd)r'>k{y) + ...|, 

where r is the Laplace sum rule variable and /, h and k are analytic functions of y. From 
this expression, one can derive the analytic form of the slope ||30|| : 



Cphysiy = 1) ~ -1 + 5,,ert + 5nP, (50) 

where at the r-stability region: 6pert — Snp — —0.04, which shows the near-cancellation 
of the non-leading corrections. Adding a generous 50% error of 0.02 for the correction 
terms, we finally deduce: 

C.,.(z/ = 1)^-1 ±0.02. (51) 

Using this result in different existing model-parametrizations, we deduce the value of the 
mixing angle: 

/I 48ds\ ^'^^ 

Vcb^i^ \ J X (37.3 ± 1.2 ±1.4) X 10-^ (52) 
where the first error comes from the data and the second one from the model dependence. 
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Let us now discuss, the effects due to the 1/M correction. It has been argued recently 
(but the situation is still controversial that the effect can lower the Isgur-Wise 

function to a value 0.89 ± 0.03 at y = 1 such that the extracted value of Vet using 
an extrapolation until this particular point will also increase by 11%. However, the data 
from different groups near this point are very inaccurate and lead to an inaccurate though 
model-independent result. Moreover, in order to see the real effect of the 1/M correction, 
one can combine this previous result at y = 1 with the sum rule estimate of the relevant 



form factor at = 0, which is about 0.53 ±0.09 With these two extremal boundary 
conditions and using the linear parametrization: 

C = Co + C'(2/-l), (53) 

we can deduce the slope: 

C' ~ -(0.72 ±0.2). (54) 

It indicates that the 1/M correction tends also to decrease which implies that, for 
larger values of y where the data are more accurate, the increases of Vcb is weaker (± 
3.7%) than the one at y = 1, which leads to the final estimate: 

-1.48psN 



X (38.8 ± 1.2 ± 1.5 ± 1.5) X 10"^ (55) 

n / 

where the new last error has been induced by the error from the slope. This result is 
more precise than the one obtained at y = 1, while the model-dependence only brings a 
relatively small error. It also shows that the value from the exclusive channels is lower 
than that from the inclusive one, which is largely affected by the large uncertainty in the 
mass definition which enters in its fifth power. Previous results for the slope and for Vcb 
are in good agreement with the new CLEO data presented in this meeting. 



7 Hybrid and Bc-mesons 

Let me conclude this talk by shortly discussing the masses of the hybrid QGQ and the 
mass and decays of the i^^-mesons. Hybrid mesons are interesting due to their exotic 
quantum numbers. Moreover, it is not clear if these states are true resonances or instead 
they only manifest as a wide continuum. The lowest cGc states appear to be a 1"*"" of 
a mass around 4.1 GeV 0. The available sum- rule analysis of the 1 state is not very 
conclusive due to the absence of stability for this channel. However, the analysis indicates 
that the spin one states are in the range 4.1-4.7 GeV. Their characteristic decays should 
occur via the rj' U(l)-like particle produced together with a ip or rjc- However, the phase- 
space suppression can be quite important for these reactions. The sum rule predicts that 
the , 0"^^ cGc states are in the range 5-5.7 GeV, i.e about 1 GeV above the spin 1. 

We have estimated the i?c-meson mass and coupling by combining the results from po- 
tential models and QSSR Potential models predict: 

Mb, = (6255 ± 20) MeV, Mb* = (6330 ± 20) MeV (56) 
MA^bcu) = (6.93 ± 0.05) GeV, Mn(bcs) = (7.00 ± 0.05) GeV 
Ms'-iccu) = (3.63 ± 0.05) GeV, M^^^bbu) = (10.21 ± 0.05) GeV 
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which are consistent with, but more precise than, the sum-rule results. The decay constant 
of the Be meson is better determined from QSSR. The average of the sum rules with the 
potential model results reads: 

(2.94 ±0.12)/,, (57) 

which leads to the leptonic decay rate into of about (3.0±0.4) x {Vcb/0-0S7y x 10^° s^^ 
We have also studied the semileptonic decay of the Be mesons and the gf^-dependence of 
the form factors. We found that in all cases, the QCD predictions increase faster than the 
usual pole dominance ones. The behaviour can be fitted with an effective pole mass of 
about 4.1-4.62 GeV instead of the 6.33 GeV one expected from a pole model. Basically, 
we also found that the each exclusive channel has almost the same rate of about 1/3 of 
the leptonic one. Detections of these particles in the next S-factory machine will serve 
as a stringent test of the results from the potential models and sum rules analysis. 



8 Conclusion 



We have shortly presented different results from QCD spectral sum rules in the heavy- 
quark sector, which are useful for further theoretical studies and complement the results 
from lattice calculations or/and heavy-quark symmetry. For the experimental point of 
view, QSSR predictions agree with available data but they also lead to some new features, 
which need to be tested in forthcoming experiments. 
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